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The structures, the electromagnetic transitions, and the beta decay strengths 
of exotic nuclei are investigated within an extended cluster model. We start 
by deriving an effective nuclear Hamiltonian within the S2 correlation opera- 
tor. Tensor forces are introduced in a perturbative expansion which includes 
up to the second order terms. Within this Hamiltonian we calculate the dis- 
tributions and the radii of A=3, 4 nuclei. For exotic nuclei characterized by 
n valence protons/neutrons we excite the structure of the closed shell nuclei 
via mixed modes formed by considering correlations operators of higher order. 
Good results have been obtained for the calculated transitions and for the beta 
decay transition probabilities. 
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1. Introduction 

The data obtained at RIA, Riken and GSI have given new life to the old 
Nuclear Physics. Nuclear structures of proton and neutron rich nuclei are 
and will be investigated giving a new insight to the fundamental observ- 
able such as the nuclear forces in the proton-proton and neutron-neutron 
components, shell closure far from stability, magnetic properties of weakly 
excited nuclear states, and many others. The theoretical analysis of these 
data requires a reliable nuclear model which can reproduce the data of 
stable nuclei and be extrapolated to predict or at least reproduce the ex- 
perimental results. The Cluster Correlation Model offers effective tools in 
this direction. We have to depart, however, from a perturbative scheme 
which is generally used to treat the two body correlation. Generalization 
to the perturbative method to include the many body correlation is in this 
paper realized within an extended, non-pcrturbative Cluster Correlation 
Model. One of the central challenges of theoretical nuclear physics is the 
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attempt to describe unknown properties of the exotic systems in terms of a 
realistic nucleon-nucleon (NN) interaction. In order to calculate matrix el- 
ements with the singular interaction (hard core) we have to define effective 
correlated Hamiltonians. 

Correlation effects in nuclei have been first introduced in nuclei by Vil- 
lars, 1 who proposed the unitary- model operator (UMO) to construct effec- 
tive operators. The method was implemented by Shakin 2 for the calculation 
of the G-matrix from hard-core interactions. Non perturbative approxima- 
tions of the UMO have been recently applied to odd nuclei in Ref. [3] and 
to even nuclei in Ref. [4]. The basic formulas of the Dynamic Correlation 
Model and of the Boson Dynamic Correlation Model (BDCM) presented 
in the above quoted papers have been obtained by separating the n-body 
correlation operator in short- and long-range components. The short-range 
component is considered up to the two body correlation while for the long 
range component the three and four body correlation operators have been 
studied. The extension of the correlation operator to high order diagrams is 
especially important in the description of exotic nuclei (open shell). In the 
short range approximation the model space of two interacting particles is 
separated in two subspaces: one which includes the shell model states and 
the other (high momentum) which is used to compute the G-matrix of the 
model. The long range component of the correlation operator has the effect 
of generating a new correlated model space (effective space) which departs 
from the originally adopted one (shell model) . The amplitudes of the model 
wave functions are calculated in terms of non linear equation of motions 
(EoM). The derived systems of commutator equations, which characterize 
the EoM, are finally linearized. Within these generalized linearization ap- 
proximations (GLA) we include in the calculation presented in the paper 
up to the ((n+l)plh) effective diagrams. The linearized terms provide, as 
explained later in the text, the additional matrix elements that convert the 
perturbative UMO expansion in an eigenvalue equation. The n-body ma- 
trix elements needed to diagonalize the resulting eigenvalue equations are 
calculated exactly via the Cluster Factorization Theory (CFT). 5 

2. The S2 correlated Hamiltonian 

In order to describe the structures and the distributions of nuclei we start 
from the following Hamiltonian: 

H = ^2(a\t\(3) alap + ^ (^ a p\vi 2 \%s) a^a^ (1) 

af3 otfi-yS 
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where V12 is the singular nucleon-nucleon two body potential. Since the two 
body states \a/3) are uncorrelated the matrix elements of v\ 2 are infinite. 
This problem can be avoided by taking matrix elements of the Hamiltonian 
between correlated states. In this paper the effect of correlation is intro- 
duced via the e tS method. In dealing with very short range correlations 
only the S 2 part of the correlation operator needs to be considered. 

Following Ref. [2] we therefore calculate an "effective Hamiltonian" by 
using only the ^2 correlation operator obtaining: 

Heff = e~ iS 'He iS ' = E Q/3 («l^)4«/3 + E^^^I^H^^ 
= Ea/jH^Ha/S + Ea/3 7 «5(*a/3bl*75)44 a ' 5a 7 

(2) 

where v[ 2 refers to the long-range part of the nucleon-nucleon force diagonal 
in the relative orbital angular momentum, after the separation: 6 

V12 = v{ 2 + v[ 2 (3) 

The separation is made in such a way that the short range part produces 
no energy shift in the pair state. 6 In doing shell model calculation with the 
Hamiltonian Eq. (2), we remark: a) only the long tail potential plays an 
essential role in the calculations of the nuclear structure i.e.: the separation 
method and the new proposed vi ow -k 7 method show a strong analogy and 
b) the v must be included as an additional re-normalization of the effective 
Hamiltonian Eq. (2). 

In Eq. (2) the ^> a( 3 is the two particle correlated wave function: 

* Q /3 = e* s *<t> a0 (4) 

In order to evaluate the effect of the tensor force on the ^l a p we calculate: 

w{r) = V^-^u(r) = V£ d ^\(nlS), J' : NL, J) (5) 

where Q is a momentum dependent projection operator, AE(ki, k 2 ) the 
energy denominator and nl the correlated two particle state in the relative 
coordinates. In Eq. (5) u(r) is generated as in Ref. 2 by a separation distance 
calculation for the central part of the force in the 3 5i state. The wave 
function obtained in this way (full line) heals to the harmonic-oscillator 
wave function (dashed line) as shown in Fig. 1. The result obtained for 
Eq. (5) calculated with the tensor force of the Yale potential 8 is given 
also in Fig. 1 left where we plot for the harmonic oscillator size parameter 
b=1.41 fm: 



*(r) = [u^Y^n^+w^Y^np)} 



(6) 



February 1, 2008 14:19 



WSPC - Proceedings Trim Size: 9in x 6in proc.TKDF 



4 




Fig. 1. Left: The u(r) and w(r) wave functions of the dcuteron, with quantum numbers 
3 Si and 3 Di, plotted as function of r; Right: Distributions of 3 H and 3 He. 



Being the admixture of the two components, circa 4%, the wave function 
Eq. (6) can be associated to the deuteron wave function. Let us use then 
the Hamiltonian Eq. (2) to calculate the structure of the A=3 nuclei. Here 
we propose to calculate the ground state of 3 H, 3 He, and 4 He within the 
EoM method which derive the eigenvalue equations by working with the 
e «s , 2 operator on the wave functions of the A=3, 4 nuclei Ref. From the di- 
agonalization of the eigenvalue equation of the three particles, we obtain an 
energy difference AE( 3 H — 3 He)=0.78 MeV and the distributions and radii 
given in Fig. 1 Right. By extending the commutator to a four particle state 
we obtain for the ground state of 4 He the binding energy of E=28.39 Mev 
and the rms radius of 1.709 fm. In dealing with complex nuclei however the 
(Si, i — 3- ■■n) correlations should also be considered. The evaluation of 
these diagrams is, due to the exponentially increasing number of terms, dif- 
ficult in a perturbation theory. We note however that one way to overcome 
this problem is to work with e^ Sl+S2+Ss ^ l- s <) operator on the Slater's 
determinant by keeping the n-body Hamiltonian uncorrelated. Via the long 
tail of the nuclear potential the Slater determinant of the "n" particle sys- 
tems are interacting with the excited Slater's determinants formed by the 
("n" particles+(mp-mh) mixed- mode excitations). The amplitudes of the 
different determinants are calculated via the EoM method. After having 
performed the diagonalization of the n-body Hamilton's operator we can 
calculate the form of the effective Hamiltonian which, by now, includes the 
complete set of the commutator equations. The method is here applied to 
6 Hc, n Bc, 14 C, 15 0, and 17 0. A detailed formulation of the model my be 
found in Ref. 9 
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3. Results 

In order to perform structure calculations for complex nuclei, we have to 
define the CMWFs base, the "single-particle energies" and to choose the 
nuclear two-body interactions. The CMWFs are defined as in Ref. 9 by in- 
cluding mixed valence modes and core-excited states. The base is then 
orthonormalizcd and, since the single particle wave functions are harmonic 
oscillators, the center-of-mass (CM) is removed. The single-particle ener- 
gies of these levels are taken from the known experimental level spectra 
of the neighboring nuclei. For the particle-particle interaction, we use the 
G-matrix obtained from Yale potential. 10 These matrix elements are evalu- 
ated by applying the e s correlation operator, truncated at the second order 
term of the expansion, to the harmonic oscillator base with size parameter 
b=1.76 fm. As elucidate in Refs. [3] and [4] the effective two-body poten- 
tial used by the DCM and the BDCM models is separated in low and high 
momentum components. Therefore, the effective model matrix elements 
calculated within the present separation method and those calculated by 
Kuo 7 in the vi ow -k approximation are pretty similar. The adopted sepa- 
ration method and the vi ow -k generate two-body matrix elements which 
are almost independent from the radial shape of the different potentials 
generally used in structure calculations. 

The particle-hole matrix elements could be calculated from the particle- 
particle matrix elements via a re-coupling transformation. In this contribu- 
tion we present application of the S„ correlated model to the charge distri- 
butions of 6 He, 11 Be, and to the electromagnetic transitions of neutron rich 
Carbon and Oxygen isotopes. The beta decay strengths from the ground 
state 14 N to the excited states of 14 C are also calculated. In Fig. 2) Left 
three distributions are given for 6 Hc: 1) the correlated charge distribution 
calculated with the full 53 operator, 2) the correlated charge distribution 
calculated with the partial S3 operator obtained by neglecting the folded 
diagrams, 3) the charge distribution calculated for two correlated protons 
in the lsi shell. The full S3 correlation operator therefore increases the 
calculated radii. In Fig. 2) Right the charge distribution for 11 Be is given. 
A charge radius of 3.12 fm has been obtained. Calculations are performed 
in a mixed S3 and S5 system. The results obtained for the Carbon and 
Oxygen isotopes are in the following presented as function of the increasing 
valence neutrons. It is worthwhile to remark that the high order correlation 
operators generate the interaction of the valence particles with the closed 
shell nucleus. The correlation model treats therefore consistently the "A" 
particles of the isotopes. 
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Fig. 2. Left:Charge distributions of °He calculated in different approximations; Right: 
Charge distribution of 11 Be. 

By using generalized linearization approximations and cluster factorization 
coefficients 5 we can perform exact calculations. In following Tables an over 
all b=1.76 fm has been used. 

In Table I, 3) we give the calculated magnetic moments and rms radii 
for one- hole and for one-particle in 16 O. The energy splitting between the 
ground- and the second (first) excited states and the electromagnetic tran- 
sitions for the two isotopes are given in Tables 2, 4). 

Tabic 1: Magnetic moment (nm) and rms (fm) of the ground state 
of ls O with J = §~;T = i 

PCM Exp. 11 

Magnetic Moment (mm) TfO .7189 

PCM Exp. 1 ^ 

rms (fm) 2.74 2.73(3) 

Table 2: Energy splitting between the ground and the second excited states 
and the corresponding electromagnetic transitions for ls O. 



Energy (McV) 


PCM 


Exp. 11 




5.41 


5.24 


Ratio 


PCM 


Exp. 11 


-B-E(-E3;| + ^i~) 


.15 


.10 


BE(M2;| + — i~) 



Tabic 3: Magnetic moment (nm) and rms (fm) of the ground state 
of 17 with J = | + ;T = \ 





PCM 


Exp. 11 


Magnetic Moment (nm) 


-1.88 


-1.89 




PCM 


Exp. 1J 


rms (fm) 


2.73 


2.72(3) 



Table 4: Energy splitting between the ground- and the first excited 
states and the E2 transition for O. 
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Energy (McV) 


DCM 


Exp. 11 




0.87 


0.89 


3 ?l '1 A 

Transition(e fm ) 


DCM 


Exp. 11 


BE(E2; \ + -> f + ) 


2.10 


2.18±0.16 



In Table 5) we give the calculated results for the energy splitting be- 
tween the ground- and the 2+ excited state and the corresponding electro- 
magnetic transition for the 14 C. The commutator equations involve S2 and 
5*3 diagrams. 

Table 5: Calculated energy splitting and BE(E2;2+ -> 0+) transition for 14 C 



Energy (McV) 


Rcf. 13 


BDCM 


Exp. 14 


A£„+,+ 




8.38 


8.32 


Transit ion (e~/m 4 ) 


Rcf. Ia 


BDCM 


Exp. 14 


BE(E2; 2 + -> 0+) 


3.38 


3.65 


3.74 ± .50 



In Table 6) preliminary results for the calculated reduced transition proba- 
bilities from the ground state of 14 N to the + , 1 + , 2 + excited states of 14 C 
are given. The calculated strengths reproduce reasonably well the experi- 
mental values. 15 

Table 6: Calculated energies of the low-lying states (MoV) + , 1+, 2+ of 14 C 
and the associated reduced transition probabilities B(GT) 
from the J=l+ T=0 ground state of 14 N. 

14 N J+T 14 C J±T Energy (McV) B(GT) 



1^0 0*1 O0 O06 

0+1 7.81 0.15 

1+1 12.17 0.12 

2+1 7.38 0.42 

2+1 8.38 0.50 

2+1 10.91 0.35 



Good results have been overall obtained for the transitions with a neu- 
tron effective charge varying between 0.1- to 0.12-e„. 

4. Conclusion and Outlook 

In this contribution we have investigated the effect of the microscopic corre- 
lation operators on the exotic structure of the Carbon and Oxygen isotopes. 
The microscopic correlation has been separated in short- and long-range 
correlations according to the definition of Shakin. The short-range correla- 
tion has been used to define the effective Hamiltonian of the model while 
the long-range correlation is used to calculate the structures and the distri- 
butions of exotic nuclei. As given in the work of Shakin, only the two-body 
short-range correlation need to be considered in order to derive the effec- 
tive Hamiltonian especially if the correlation is of very short range. For the 
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long range correlation operator the three body component is important and 
should not be neglected. Within the three body correlation operator, one 
introduces in the theory a three body interaction which compensates for 
the use of the genuine three body interaction of the no-core shell model. 
Within the S2 effective Hamiltonian, good results have been obtained for 
the ground state energies and the distributions of 3 H, 3 He, and He. The 
higher order correlation operators S = 3 • • • n have been used to calculate 
the structure and the electromagnetic transitions of ground and first ex- 
cited states for the isotopes of Carbon and Oxygen. By using generalized 
linearization approximations and cluster factorization coefficients we can 
perform expedite and exact calculations. 
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